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PREFACE 


Part  of  the  RAND  research  program  consists  of  basic  sup¬ 
porting  studies  in  mathematics.  The  present  paper  continues  the 
work  of  the  authors  of  invariant— imbedding  techniques  for  solving 
problems  in  radiative  transport  theory. 


SUMMARY 


In  the  present  paper  it  is  shown  how  the  invariant— imbedding 
technique  is  used  for  the  derivation  of  integral  equations  governing 
the  reflection  and  the  transmission  coefficients  of  radiation  in  a 
one- dimensional  medium  with  a  moving  boundary.  In  Part  A,  when 
formulating  the  source  function,  account  is  taken  by  the  relaxation 
of  photon  emission  under  the  assumption  that  during  an  elementary 
act  of  scattering  a  photon  remains  in  a  state  of  absorption  for  a 
certain  duration  of  temporal  capture.  In  Part  B,  the  Schuster 
problem  for  a  moving  one-dimensional  medium  is  exactly  treated. 
In  other  words,  when  formulating  the  transfer  equations,  the 
influence  of  the  Doppler  effect  caused  by  the  macroscopic  dif¬ 
ferential  motions  is  considered  as  a  principal  agency  of  the  change 
of  frequency  on  scattering. 
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INVARIANT  IMBEDDING  AND  SCATTERING  OF  LIGHT  IN  A 


ONE-DIMENSIONAL  MEDIUM  WITH  A  MOVING  BOUNDARY 

1.  INTRODUCTION 

In  recent  years,  for  analytical  and  numerical  advantages,  the 
principle— ofHLnvariance  method  originally  due  to  Ambarzumian  [l] 
and  Chandrasekhar  [6]  has  been  developed  by  several  authors 
(Bellman,  Kalaba,  and  Wing  [2]  and  [3]  ,  Wing  £  11] ) .  The 
invariant— imbedding  technique  thus  proposed  was  successfully 
applied  not  only  to  various  kinds  of  stationary  transport  problems 
for  photons  and  neutrons,  bu+  also  to  time— dependent  transfer 
problems.  Many  previous  attempts  at  solving  nonsteady  transfer 
problems  will  be  found  in  the  list  of  references  of  the  preceding 
papers  (see  also  Bellman,  Kalaba,  and  Ueno  [4]  and  [5]  ,  Kaplan 
[  8  ]  ,  and  Wing  [  1 0  ] ) . 

Recently,  in  the  study  of  astrophysics,  emphasis  has  been 
put  on  the  elucidation  of  the  photon  diffusion  process  in  a  dif¬ 
ferentially  moving  atmosphere,  whereas  the  exact  treatment  seems 
very  difficult  on  account  of  the  untractability  of  the  transfer 
equation  due  to  the  Doppler  shift  in  frequency.  Such  problems  are 
encountered  in  the  interpretation  of  spectral  lines  in  novae, 
planetary  nebulae,  and  the  solar  prominences  and  corona.  This 
state  of  affairs  is  also  found  in  the  analysis  of  several  physical 
experiments . 
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By  reducing  the  transfer  problem  to  a  novel  boundary-value 
problem  for  hyperbolic  differential  equations,  Chandrasekhar  [7] 
has  given  approximately  the  solutions  illustrating  the  effects  that 
may  be  expected  in  the  rectangular  contour  of  the  scattering  coef¬ 
ficient.  Later,  combining  the  Doppler  effect  due  to  the  thermal 
motions  of  atoms  and  the  macroscopic  differential  motions  in  finite 
one—  and  three-dimensional  media  and  assuming  an  arbitrary 
contour  of  the  scattering  coefficient,  Sobolev  [9]  has  considered 
the  approximate  effect  of  the  complete  redistribu'ion  of  radiation  in 
frequency  on  the  radiation  pressure  due  to  the  Lyman— a. 

In  a  preceding  paper,  the  invariant— imbedding  technique  was 
applied  to  the  study  of  time-dependent  scattering  of  light  in  a  one- 
dimensional  medium  of  a  finite  optical  thickness.  In  this  paper, 
extending  the  imbedding  procedure  to  a  one— dimensional  moving 
medium,  we  shall  consider  a  couple  of  nonsteady  scattering  pro¬ 
cesses  as  follows:  In  Part  A  the  left-hand  boundary  of  the  medium 
moves  as  a  prescribed  function  of  time,  and  the  source  function  is 
formulated  under  consideration  of  the  relaxation  of  the  photon 
emission.  In  Part  B,  we  treat  the  Shuster  problem  for  a  dif¬ 
ferentially  moving  atmosphere  such  that,  assuming  that  the  only 
optical  effect  of  large-scale  motions  is  a  change  of  frequency  on 
scattering,  the  interaction  of  radiation  fields  with  different  fre¬ 
quencies  is  taken  into  account.  In  a  later  paper,  extending  this 


I 
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procedure  and  allowing  for  the  noncoherent  scattering,  we  shall 
consider  the  formation  of  spectral  lines  in  a  differentially  moving 
atmosphere  in  rod  and  slab  geometries. 

A.  RELAXATION  OF  PHOTON  EMISSION 
2.  THE  EQUATION  OF  TRANSFER 

Consider  a  one-dimensional  medium  of  finite  geometrical 
thickness  z^  -  zQ  >  0  with  a  moving  boundary  whose  optical 
properties  vary  with  its  position  and  with  time.  Let  a  pencil  of 
radiation  of  Dirac  delta-function  time-dependent  intensity 
XQ(t)  =  F  6(t  -s)  (t  >  s  >0)  be  incident  on  the  fixed  boundary 
where  F  is  a  constant  and  6  is  the  Dirac  delta-function  (see 
Fig.  1).  At  time  t,  the  position  of  the  left-hand  boundary  is 
zQ  =  zQ(t),  whereas  at  t  =  0  it  is  equal  to  zero;  further,  |dzQ/dt| 
is  less  than  the  speed  of  light  in  the  medium. 


zo  =  0 


t  =  0 


*o"2o(8“T)  t  =  s‘t 


z.  +  A 

1  Incident 
radiation 

— I  < - 


z^+  A 


Transmitted 
radiation 
« - 


I  (z,  t)  I  (z,t) 

l - — y 


zo=zo<° 


t  >  s 


Reflected 

radiation 

- > 


z^  +  A 


Fig.  1  The  reflection  and  the  transmission  of  radiation  in  a 
one— dimensional  medium  with  a  moving  boundary 


i 
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For  the  sake  of  simplicity,  the  scattering  of  light  in  either 
direction  is  assumed  to  be  equally  probable. 

Let  Ij(z,  t)  denote  the  specific  intensity  of  radiation  at  depth 
z,  at  time  t,  directed  toward  the  boundary  7.  and  similarly  let 
^  denote  that  directed  toward  the  boundary  Zg(t). 

The  equation  of  transfer  appropriate  to  the  present  case  is 
written  in  the  form 

9I1  1  9I1 

(2>1)  -aT  +  c“ar  =  -i(z’  tni{z>  t>  +  fi(z>  t), 
ai2  1  ai2 

^2‘  3z  +  c  3z  ^  + 

where  c  is  the  speed  of  light,  1  is  the  extinction  coefficient,  and 
B(z,  t)  is  the  source  function. 

In  a  manner  similar  to  that  used  in  a  preceding  paper  [5]  , 
the  source  function  takes  the  form 

(2.3)  B(z,  t)  =  ~  a(z,  t)  J  ft  (z,  t')  +  I  (z,  t')}  e"^ 

*s 

where  a  is  the  scattering  coefficient  and  t  is  the  duration  of 
temporal  capture.  Equation  (2.  3)  represents  the  relaxation  of 
emission  of  absorbed  energy. 

Equations  (2.  1)  and  (2.2)  should  be  solved  subject  to  the 


boundary  and  initial  conditions 
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(2.4)  Ii(zQ(t),  t)  =  0,  I2(Zl,  t)  =  F6(t  —  s), 

(2.  5)  I^z,  s)  =  0,  I2(z,  s)  =  0  (zQ  <  z  <  Zj,  0  <  s  <  t), 

where  I^(z^,  0  and  I2(zQ,  t)  are  respectively  the  reflected  and  the 
transmitted  intensities. 

3.  THE  REFLECTED  INTENSITY 

Let  the  coefficient  of  reflection  be  denoted  by  R(z^;  t,  s). 

Then  the  reflected  intensity  is  expressed  in  terms  of  the  reflection 
coefficient  R  as  follows: 

(3.1)  I1(z1,  t)  J  R(zx;  t,  t')I2(z1,  t')dt', 

s 

where  I,,^,  ^  is  given  by  (2. 4). 

We  shall  inquire  into  an  integral  equation  for  R(z^;  t,  s) 
with  the  aid  of  the  invariant— imbedding  technique. 

When  we  add  a  layer  of  infinitesimal  thickness  a  to  the 
boundary  z^,  the  imbedding  procedure  of  the  scattering  in  the  added 
layer  yields 

(3.2)  I^Zj+A,  t  +  *)  =  I*(Zl,  t)+  {-f(zr  t)  I1(z,  t)+B(zr  t)}A+(P(A), 

where  0(A)  is  of  the  order  of  magnitude  of  the  infinitesimal  A2. 
Equation  (3.  2)  shows  that  the  reflected  intensity  at  the  boundary 
z  +  A  at  time  t  +  — ,  due  to  the  incident  intensity  I  (z  +  A,  t1), 

A  C  Ct  1 
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is  produced  by  the  intensity  at  the  boundary  z at  time  t,  due  to 
a  modified  incident  radiation  source  ^2(z  ,  addition  to  the 

intensity  produced  by  the  interaction  of  the  radiation  field  in 
(zr  Zj  +  A). 

The  term  I  on  the  left-hand  side  of  (3.  2)  is  given  by 

t+±- 

r  c  . 

(3.3)  I^Zj  +  A,  t  +  —)  =  J  R(z1  +  A;t  +— ,  x)I2(Zl  +  A,  x)dx, 

s  —  — 
c 

where  I^z^  +  A,  t)  =  F  6(t  —  s  +  ^). 

On  the  other  hand,  recalling  (2.  2),  we  obtain  the  modified 
incident  source  given  by 


(3.4)  I2(Zl,  t)  =  I2(Zl  +  A,  t  -*)  -i( zv  t)I2(Zl,  t)i 

t  — (t  —  x)/t 

+  ^-ct(z1,  t)  Ajf  {I1(z1,  x)  +  I2(z1,  x)}  e 


1  dx 


t$(A). 


Thus,  the  intensity  of  radiation  due  to  the  modified  source 


l2  iS 


*  A  A 

(3.5)  I^z^  t)  =  /  R(z1;  t,  x)  I2(z1,  x)dx 

t  S  4 

=  J  R(zx;  t,  x)I2(z1+A,  x  —  ~)dx  —  A  /  Rtz^  t,  x)l(z1,  a)  1^,  x)dx 


i  r1  rx 

j  A  J  R(z:;  t,  x)  V(zv  x)dx  j  I^Zj,  y> 


-(x-y)/t 


ldy 

t. 
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A J'  R(z1;  t,  x)a(z1, 


x 

12{ZV  y)  6 


-(x  —  y)/t 


1 


2 


On  substituting  (3.  3)  and  (3.5)  into  (3.2)  and  putting  I2(z  ,  ^ 


IQ(t),  we  get 


(3.6)  R(z1  +  A;  t  +  ^,  s  “■  “)  =  Rtz^  t,  s)  -  i( z^,  s)  R(z;  t,  s)A 


-<!i-s)/,l  dx 


+  f  A  /'r(z1;  t,  x)d{ z1,  x)dx  J  R(z,;  y,  x) 

1  r* 

+  2AJ  R(zi:  t>  xMz^  x) 

+  A|j-i(z1>  t)  R(Zjj  t,  s)+|-a(z1,  t)^ R(Zjj  x,  s) 

1  +SXM. 


{x-y)/tl  ^ 


— (t  —  x)/t 


1  dx 
t. 


-  (t  -  s)  / 1 . 


+  i<7(z,,  Of 


Then,  in  the  limit  as  A  -*  0,  (3.6)  becomes 


<3'7)  +  t)  +  «v  S»R 

—  (t  —  sJ/tj  t 

=  2  a(zl»  f  +  i  a*zi’  ^ J'  R^Z1;  x>  s) 


~(t  -x)/t 


1  dx 
t. 


1  / 

+  R(z^  t,  x)  a(zx,  x) 


-  (x  -  x)/t 


1  dx 
t. 


i  rx  -(x-y)/t1  . 

J y  R(z^  t,  x)  a(Zj,  x)dx y  R(zx;  y,  s)  e 


Equation  (3.  7)  is  the  requisite  integral  equation  for  the  reflection 
coefficient.  The  boundary  and  initial  conditions  imposed  on  R  are 
given  by 
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(3. 

8) 

R(z1j  t, 

s)  =  0 

for 

zo(s)  =  Z1 

(3. 

9) 

R(z1;  t, 

s)  =  0 

for 

s  >  t  >  0. 

When  the  optical  properties  of  the  medium  are  steady  and 
homogeneous  with  respect  to  time  and  depth,  and  we  put 


(3. 10) 


T  =  iz, 


7 


where  t  is  the  optical  depth,  t ^  is  the  mean  free  time,  and  25 
is  the  albedo  for  single  scattering,  (3.  7)  reduces  to 


(3.„,  jS- (Tl;  t  -  ,>  ♦  ,2  |5_l2|5t2R 


-  (t  -  s)/t 


1 

+  2 


rv  -(x-s)/t 

'■J  R(ti;  t  —  x)  e  1-2* 


■Ifc 

s 

rl  rx  “(x-y)/t  -| 

+J  R{Ti‘ 1  ~ x)dx  J  R(iv  y_s) e  1  "fj  • 


4.  THE  TRANSMITTED  INTENSITY 

Consider  the  case  of  illumination  of  the  fixed  boundary  z  . 
Let  the  coefficient  of  transmission  be  denoted  by  T(zJ;  zQ(t),  t; 
Zq(s),  s).  Then  the  transmitted  intensity  at  the  boundary  zQ  at 
time  t  is  expressed  in  the  form 

rt* 

(4.1)  I2(z0(t),  t)  =J  T(z1;  zQ(t),  t;  zQ(x),  x)  x)dx, 

s 

where  I2(z1,  t)  is  given  by  (2. 4)  and  t*  =  t  -  -  zQ  (t))/c. 


I 
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For  convenience,  we  shall  write 
(4.2)  T(z1;  zQ(t),  t;  zQ(s),  s)  =  T(z1;  t,  s). 

In  what  follows,  we  shall  seek  an  integral  equation  for  the 
T— function.  On  adding  a  layer  of  infinitesimal  thickness  A  to  the 
boundary  z^,  we  get  (3.4).  On  applying  the  imbedding  argument 
to  the  transmitted  intensity,  we  obtain 


(4.3)  I^Zgtt),  t)  =l2(zQ(t),  t). 


On  the  other  hand,  we  have 


t*  —  &lc 

(4.4)  I2(zQ(t),  t)  *  J  T(z1+A;t,  x)  I^Zj  +  A,  x)dx, 


and 


s  —  ■ 


,t* 


(4.5)  I2(zQ(t),  t)  =  P  T(z1;  t,  x)/I2(z1,  x)dx. 

Substitution  of  (3.  4)  into  (4.  5)  yields 


(4.6)  I2(zQ(t),  t)/F  =  T(zx;  t,  s)  -Kz^  s)  Tfz^  t,  s)A 

-•*  -(x-s),tldx 

*1 


t* 

2  A  J  T(z1;  t,  x)a(z1,  x)  e 

^  c* 


.t* 


i  r  rx  -(x-y)/t -  . 

+  j  A  J  T(z1;  t,  x)  o(z1,  x)dx  J  R(z1;  y,  s)  e  +(9(A). 


On  combining  (2.4)  and  (4.3)— (4.5),  we  get 


/ 
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(4.  7)  T(Zl  +  A;  t,  s  -=)  =  Tfz^  t,  s)  -ifz^  s)  t,  s)A 

li 


t* 

+  j  A J  T(Zl;  t,  s)  a(Zl,  x)  e 


±a/  T(Zl;  t,  x)  a(z1?  x)dx y,  s) 


-(x-y)/t1 


‘l 


(9(A). 


Letting  A  -*  0  in  (4.6)  gives 


<4-8>  + 

Jt* 


-  (x  -  s) /  t 


1  dx 

t. 


~  \  t>  xjfftej,  x) 

i  rt#  rx  -(x-yj/tj 

+  -J  T(Zl;  t,  x)  a(zi;  x)dx J  Rtz^  y,  s)  e  -p 

Equation  (4.  8)  is  the  required  integral  equation  governing  the 
transmission  coefficient,  where  T(z^;  t,  s)  =  0  for 


t  —  s  < 


zi  _zo(t) 


Let  the  reduced  and  the  diffuse  components  of  the  transmission 
coefficient  be  denoted  by  T°(z1;  t,  s)  and  T*^;  t,  s),  respectively, 
where 


(4.9) 


r  ri  z  —  z  -[  z  —  z  (t) 

T°(z1;  t,  s)  =  exp  j—  J  i(z,  s  + - - —  dzj  6(t - - - s) 


zo(t) 
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Then  the  insertion  of  T  =  T°  +  T*  into  (4.  8)  provides  us  with 
an  integral  equation  for  the  T*— function. 

In  the  steady  and  homogeneous  case  with  respect  to  time  and 
depth,  by  (3.10),  (4.8)  becomes 


(4.10)  £L-(Tl;.-s)-,2g  +  T 


■f  [f  T(Ti ;t_x) 


-  (x  -  s) / t 


1  dx 
t. 


rx 

+  /  Tfr^t-xjdx/  Rfi^y-s) 


-(x-y)/t1 


dyl 

V’ 


where  t*  =  t  —  (t^  —  t^)^. 

B.  CHANGE  OF  FREQUENCY  ON  SCATTERING  DUE 
TO  DOPPLER  EFFECT 

5.  THE  EQUATION  OF  TRANSFER 

Consider  a  one-dimensional,  differentially  expanding  medium 

of  finite  geometrical  thickness  z^  —  z^>  0  whose  velocity  v(z)  >  0 

at  z  is  a  linear  function  of  the  depth,  z^  <  z  <  z^.  ^et  a^v>  z ) 

denote  the  monochromatic  volume  scattering  coefficient  at  z  for 

the  frequency  v  as  judged  by  a  stationary  observer.  The  observer 

may  be  considered  to  be  at  rest  with  respect  to  the  position  at  z  =  0. 

For  simplicity,  the  quantity  a  is  assumed  to  be  p(z)  (7*(v),  where 

p(z)  is  the  mass  density  and  <7*(v)  is  the  mass  scattering  coefficient. 

Furthermore,  we  assume  the  characteristic  feature  of  a*,  that  the 

only  optical  effect  of  the  large-scale  motions  v(z)  is  a  change  of 

frequency  on  scattering. 
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Let  I^v,  z)  denote  the  specific  intensity  of  radiation  at  z 
directed  toward  the  boundary  z  and  of  frequency  v  as  judged  by 
the  stationary  observer,  and  further  let  I2(v,  z)  denote  the  specific 
intensity  directed  toward  the  boundary  z^  (see  Fig.  2). 


Incident 
radiation  ^ 


Reflected  z0  ^  Z0 
radiation 


I  (v,  z)  lAv,  z)  Transmitted 

4: — — -  — — - — >  radiation 


Fig.  2.  The  reflected  and  transmitted  intensities  of  radiation 
by  a  differentially  moving  one-dimensional  medium 


Consider  a  pencil  of  radiation  of  constant  intensity  F  in  all 
frequencies  to  be  incident  on  z^.  When  we  take  into  account  the 
Doppler  shift  in  frequency,  -^v(z)(l  +  l),  caused  by  the  large- 
scale  motions,  the  equations  of  transfer  appropriate  to  the  present 
case  take  the  form  (cf.  Chandrasekhar  [7]) 


(5.1) 


(5.2)  - 


9I1(v,  z) 
9z 

9I2(v,  z) 
9z 


-|ct(v(1-^),  z)  [Ij^v,  z)-I2(v(l— 2^),  z)], 
-■|  CT(v(l  +  ^ ) ,  z)  [I2(v,  z)  -IjMl  +  2^),  z)]  , 


where  c  is  the  speed  of  light. 

In  what  follows,  in  a  manner  similar  to  that  used  by 
Chandrasekhar  [7]  ,  we  shall  make  some  simplifications  of  the 
transfer  equations.  In  the  evaluation  of  the  Doppler  shift  in  frequency, 
v(l  —  ^)  is  reasonably  approximated  by  v  —  where  is  the 
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frequency  of  the  line  center. 

Furthermore,  in  place  of  considering  the  equations  of  transfer 
for  some  fixed  frequency  v,  we  shall  somewhat  modify  them  in  such 
a  way  that  v  in  (5. 1)  and  v  in  (5.  2)  are  respectively  replaced 
by  Vj  and  v 2,  where 

(5.3)  Vj  =  v  +  vQ  v2  =  V-V0c- 

Allowing  for  (5.3),  (5.1)  and  (5.  2)  are  rewritten  as  follows; 

31 

(5.4)  ~  (vx,  z)  =-|  zHVV  z)  _I2(v2’  z)1  » 

01 

(5.5)  - —  (v2,  z)  *  -|cr(v,  z)tI2(v2'  ->-VV  Z)^  ' 

Writing 


(5.6)  I^v.,  z)  =  «t(v,  z)  (i  =  1,  2), 


we  have 


81  9^,  _  9^ 

(5-7)  Z)  =  ~+Q“ 


for  i  =  { 2  , 


where  Q  is  a  constant  given  by 


(5.8)  Q=  — 


0  dv(z) 


By  (5.  7),  (5.4)  and  (5.  5)  become 
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(5.9) 

~J  <7(v,  z)  [ -  *2]  , 

3z  ^  3v 

(5.  10) 

3 

3z  +  Q  3v  = 

-|ct(v,  z)  [^  -  *2]  . 

Equations  (5.9)  and  (5.  10)  should  be  solved  subject  to  the 
boundary  conditions 

(5.11)  ^(v,  Z(J)»F,  *2(v,  Z]L)  =  0, 

where  F  is  a  constant  for  all  frequencies. 

The  quantities  I2(v,  zg)  and  Ij(v,  z^)  represent  respectively 
the  reflected  and  the  transmitted  intensities. 

6.  THE  REFLECTION  COEFFICIENT 

Let  the  coefficient  of  reflection  be  denoted  by  R(Zg,  z^;  v). 
Then  the  reflected  intensity  ^(v,  zq*  expressed  in  terms  of 
R  as  follows: 

(6.  1)  *2(v,  zQ)  =  F  R(z0,  zx;  v). 

In  a  manner  similar  to  that  used  in  Sec.  3,  we  shall  seek  an 
integral  equation  for  the  R— function  by  means  of  the  invariant- 
imbedding  technique. 

By  addition  of  a  layer  of  infinitesimal  thickness  A  to  the 
boundary  z^,  and  use  of  (5.10),  the  imbedding  procedure  gives 
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(6.2)  *2(v-QA,  zQ  -  A)  3  **(v,  zq)+|-  ff(v,  zq)A(V£v,  zq)  -  z^  +$A). 

By  (5.9),  the  modified  intensity  at  z  in  an  added  layer  is 
provided  by 

(6.3)  zQ)  =¥j(v  +QA,  zQ— A)  — |<T(v,  z^a[^(v,  zq)  z^]  +  (Ra). 

On  making  use  of  (6.  3),  we  see  that  the  reflected  intensity  of 
radiation  due  to  the  modified  source  ^  becomes 

(6.4)  ¥*(v,  zQ)  =  F  [R(zQ,  z^  v)  - 1  CT(v,  zQ)  R(zQ,  z^  v)A 

+  |  CT(v,  zq)  R2(zq,  ZjJ  v>\]  +(P(A). 

The  combination  of  (6.  1),  (6.2),  and  (6.4)  provides  us  with 

(6.5)  R(zq  —  A,  ZjJ  v  — QA)  =  R(zq,  z1»  v) 

+  |  a(v,  zq)A  — <t(v,  zq)R(z0,  ZjJ  v)  A  +  |a(v,  zQ)R2(z0,  Zjj  v)A  +^A). 
Then,  letting  A  -  0,  we  have 

(6  6)  +q|5  =-|a(v,  z0)+C7(v,  zq)R 

-|a(v,  zq)R2. 

Equation  (6.6)  is  the  required  integral  equation  for  the 


reflection  coefficient  R. 
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7.  THE  TRANSMISSION  COEFFICIENT 

Consider  the  transmitted  intensity  of  radiation  emergent  from 
the  boundary  z^  for  the  illumination  of  the  boundary  z^. 

Let  the  coeff’"  lent  of  transmission  be  denoted  by  T(z^,  z^;  v). 
The  transmitted  intensity  at  z^  is  expressed  in  the  form 

(7.1)  zx)  =  F  T(Z(),  zx;  v). 

On  adding  a  layer  of  infinitesimal  thickness  A  to  the  boundary 
Zg,  we  get  (6.  3).  By  imbedding  the  photon  diffusion  process  in  the 
added  A— layer,  we  obtain 

(7.2)  *J(v,  Zj)  =  **(v,  zx). 

On  the  other  hand,  we  have 

(7.3)  *°(v,  Zj)  =  F  T(zQ  -  A,  Zy  v  +QA), 
and 

(7.4)  tf*(v,  Zj)  =  T(Zq,  Zy  v)  ^x(v,  zQ). 

Then,  combining  (6.  3)  and  (7.  3)— (7.4),  we  obtain  (7.2)  in  the 

form 


T(zq  -  A,  Zy  v  +  QA)  =  T(zq,  z^;  v)  [  1  — 

-|a(v,  Zq)  A  {l  —  R(zq,  Zy  v)}]  +(9(A). 


(7.5) 
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In  the  limit  as  A  -*•  0,  we  have 

ftT1  ST*  1  1 

<’•«  =  2  <7^V'  Z0>T-2<7(v'  Z0)TR- 

Equation  (7.  6)  is  the  requisite  integral  equation  governing  the 
transmission  coefficient  T. 


/ 
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